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Abstract. In a recent paper, Yu. A. Brychkov derived a series of identities for 
multiples sums of special functions, using generating functions. Among these 
identities, a particularly interesting one involves multiples sums of Bessel I v 
functions with half-integer indices. We derive here some equivalent identities 
that involve different kinds of Bessel polynomials, sing a probabilistic approach 
based on the properties of the Generalized Inverse Gaussian probability den- 
sity. 



1. Introduction 

In a recent paper [T], using generating functions, Yu. A. Brychkov derived the 
following identity 
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where Ik (z) is the modified Bessel function of the first kind. Our aim here is to 
show that this formula is equivalent to several simple identities involving various 
kinds of Bessel polynomials, and to provide some extended versions of it. Our 
approach involves a probabilistic interpretation of these identities which does not 
require the computation of any generating function. 
Let us first remark that, from [3l p. 675], 

J_ fc _i (z)-J fc+ i (z) = l(~l) k K kH (z) 

where K k (z) is the modified Bessel function of the third kind or Macdonald func- 
tion. For k integer, this function is related to the Bessel polynomial q k (z) of degree 
k as follows 
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where the Bessel polynomial q k (z) is defined by 
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First examples of these Bessel polynomials are 



<7o {z) = 1, qi (z) = l + z, q 2 (z) = 1 + z + —. 
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These polynomials satisfy the normalization constraint 

g„(0) = i. 

Replacing the Bessel functions Ik (z) by their expression in terms of Bessel poly- 
nomials yields, after some elementary algebra, the equivalent and more compact 
version of identity ((T|) 
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Moreover, replacing the Bessel functions Ik (z) by their expression in terms of Bessel 
Kk (z) functions yields the equivalent identity 
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2. A PROBABILISTIC APPROACH TO IDENTITY © 

We show here that identity ([3]) can be interpreted in a probabilistic setting; this 
interpretation relies on properties of the generalized inverse Gaussian distribution, 
defined as [2J 

(±\^ 

(4) / (a;; i/,, X , A) = ^^^- x^ 1 exp ("fa;" 1 - -*) , x, iff, X > 0, A > -1. 

In order to explicit this probabilistic interpretation, we need the following pre- 
liminary results. This first lemma can be found in [2]. 

Lemma 2.1. |2J Let X_i Zl and X_i Z2 two independent random variables with 

generalized inverse distribution as in with parameters tp — 1, x — z 2 an d 
lambda = — |, then 

( 5 ) X -i, Zl + X -±,z 2 ~ X -±,z 1+ z 2 - 

where the sign ~ denotes identity in distribution. Moreover, 

(6) X_ hz +X h0 ~X hz . 

From these two results, we can deduce the following identity, that can also be 
found in [2] 

Lemma 2.2. Let Xi Zi and Xi Z2 two independent random variables with general- 
ized inverse Gaussian distribution as in with parameters = 1, x — z2 i A = 5, 
t/ien t/ie identity in distribution 

( 7 ) +JC |^2 ~^|,o+A^i, zl+22 
holds, where Xi Zl+Z2 is independent of Xi . 

Proof. We remark from ([S]) that 
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where the random variables on the right-hand side are mutually independent. By 
the stability property (O, we obtain 

Xi )Zl +Xi ;Z2 - X_i }Zl+Z2 +Xi fi + Xi fi 

where we have used again the property ([6]), hence the result . □ 

We note that the probability density of Xi is chi-square density with one 
degree of freedom, or equivalently a Gamma density with scale parameter 2 and 
shape parameter 1/2: 
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The link between the Generalized Inverse Gaussian random variables and the 
Bessel polynomials is characterized as follows. 

Lemma 2.3. Given v G R. the v—th moment of Xi z is 

(8) EXl = W^exp (z) z v+ *K v+ i (z) . 

When v — n is an integer, this formula simplifies to 

1 (2nV 

0) ^Xl z = - [ -± qn {z) 

where q n (z) is the Bessel polynomial of degree n. 

Proof. The v— th moment is easily computed as the integral 

x v f (x; 1, z 2 , dx = y|exp (z) z v+ ^K p+h (z) 

where Kk is the Bessel function of the second kind of order k. In the case where 
v — n is an integer, using the expression ([2]) of the Bessel function K n+ i in terms 
of the Bessel polynomial q n , we obtain 

2™ / 1\ 

with ^rrH2) — 2^"^T^' hence the result. □ 
We now have the necessary tools to prove the following extension of identity ([3]) . 
Theorem 2.4. // {zj, 1 < i < m; m > 2} are complex numbers then, with z = 

y" m z- 

( io ) e n( vfc j^^)=2:( vfc j — — »(*)■ 

The special case where all zi are equal reads Brychkov's identity (0). 
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Proof. We consider first the case to = 2. As a consequence of Lemma 12.21 the 
moments of Xi Zi + Xi Z2 are the same as the moments of Xi + Xi^ Zi+Z2 . Using 
the binomial expansion, the moment of order n of Xi Zi + Xi Z2 yields 

E(x hzi +X h S = J2(l)E(x hzi ) k E(x l2 
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The same approach applied to Xi + X± Zl + Z2 gives 
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which yields the result. The extension to the case of an arbitrary integer value 
to > 2 is left to the reader, using the following elementary extension of identity ([7]): 

Lemma 2.5. // {z^, 1 < i < m; m > 2} are real positive numbers and Xi are in- 
dependent random variables, then the following identity in distribution holds 

m 

J2 X h^- Xi ^^+X^i^. 

i=\ 

□ 

We remark that X m -i n is distributed as a Gamma random variable with scale 

2 > u 

parameter 2 and shape parameter m ~ 1 (or equivalently a chi random variable with 
to — 1 degrees of freedom) . 

3. LINKS TO LAGUERRE POLYNOMIALS 

We note that the Laguerre polynomials L„ (z) are related to the Bessel poly- 
nomials q n (z) as 

q n (z)= { -^Li- 2n -V (2z); 

as a consequence, an equivalent statement of identity (|10p in terms of Laguerre 
polynomials reads, with z = 53fc=i z ii 

E II L i~ 2ki ^ (*) = E 2 -rtkr ( - 4)n " L * 2k ~ 1] (z) • 

fclH Vk m =ni=l fc=0 
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The special case m = 2 and Z\ = — Z2 = z of this identity can be found in [H 4.4.2.9] 
as 

f:4- 2 ^ 1) w4-r 2fc - 13 (--) = (-4r- 

k=0 

4. LINKS TO OTHER BESSEL POLYNOMIALS 

In this section, we show that two other families of Bessel polynomials can be 
interpreted as moments in the same way as identity (|9j). We deduce from these 
representations some identities equivalent to (|10[) but that have a more simple 
form. 

4.1. Bessel 8 n polynomials. In his textbook [5], Grosswald considers the follow- 
ing Bessel polynomials, called reverse Bessel polynomials 

(11) e n (z) = ^2- n q n (z) . 
For example, 

O (z) = 1, 0i (z) = 1 + z, 2 (2) = 3 + 3z + z 2 . 
From identity ©, we deduce 

e n {z) = ex? z . 

Moreover, the identity (|10p reads in terms of these polynomials 

\p TT 0fcj (ft) _ Z I 2 ) n -k Ok (Z) 
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with z = Zi. The case m = 2 appears in [6l eqn. (5.4)]. 

4.2. Bessel /„ polynomials. The third family of Bessel polynomials /„ is defined 
by L. Carlitz [6] as 

(12) /„ (z) = z0 n _i (z) ; n > 1, 
and /o (z) = 1. First examples are 

fl (z) = z; / 2 (z) - z + z 2 ; / 3 (z) - 3z + 3z 2 + z 3 . 

Note that as in the case of the reverse Bessel polynomials, these polynomials have 
their highest degree coefficient equal to 1. 

Since it can be easily checked from (Q} that the Generalized Inverse Gaussian 
density satisfies the functional equation 

xf ( x; l,z 2 ,-i ) = zf (x; l,z ? 



2 J J V 2, 
we deduce 

zEX?- 1 = EX n 1 

so that 

/„ (z) = EXI\ z , Vn > 1. 

Since moreover /o (z) = 1, this representation holds in fact Vn > 0. By the stability 
property ([2]), we deduce that the polynomials /„ satisfy the multinomial property 
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with z — X)i=i z i- The case m = 2 is given in 6, eqn. (2.7)]. 

Moreover, since, by (fT2|) . (fTTj) and ([2]), the polynomials /„ are related to the 
Bcssel functions K n as 



fn (*) = 




we deduce an equivalent version of ()13[) 

N : nv*-* w -(i) *V*-»<«> 

fciH Hs m =tii=l v 7 

with 2 = 5Z i=1 2», which can be found in 3, 5.18.1.3]. 



5. Conclusion 

We have shown that an identity introduced by Brychkov using generating func- 
tions can be interpreted as a multiplication identity for several types of Bessel 
polynomials. Moreover, we have given a probabilistic background to this iden- 
tity and exhibited its relationship with the Generalized Inverse Gaussian density. 
Other tools, such that generating functions, can certainly replace this probabilistic 
approach; however, we found that it is particularly convenient in this context. 

As a final illustration of this efficiency of this tool, we derive a quick proof of a 
famous Turan-type inequality for Bessel K v functions: in [5], Ismail and Muldoon 
proved that the function 

K v {z) 

is increasing, implying that the function 

v i-> K v (z) 

is log-convex. As a consequence, the following Turan-type inequality - as named 
by Karlin and Szego - for Bessel functions holds: for - and % > — i, and with 

p q 

K± + y q (z) < K x (z)*K y (z)« . 

The case p = q = 2 of this inequality is also derived in [5] from the following 
generalization of the Cauchy-Schwarz inequality 

" g{t)r(t)dt x j\(t)P(t)dt> y\{t)f^(t)dtj 

The probabilistic approach based on the moment representation © allows to 
derive a straightforward proof of this result that does not require such a refinement: 
choosing a random variable Xi z and applying the standard Holder inequality 

E (ZiZ 2 ) < [EZf]p \EZ$\* 
with - + - = 1 and Z\ = X V J 2 ^ 4 and Zo — X^ 2 1 ^ 4 , we deduce the result. 
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